Self-pulsing of electron transmission by a transversal magnetic field 
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The distribution of scattering delay times is analyzed for classical electrons which are transmitted 
through a finite waveguide. For non-zero magnetic field the distribution shows a regular pattern of 
maxima (logarithmic singularities). Although their location follows from a simple commensurability 
condition, there is no straightforward geometric explanation of this self-pulsing effect. Rather it 
can be understood as a time-dependent analog of transverse magnetic focusing, in terms of the 
stationary points of the delay time. We also discuss the possibility of singularities in the delay-time 
distribution for generic 2D scattering systems. 
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Ballistic semiconductor nanostructures offer fascinat- 
ing opportunities to observe effects from the theory of 
classical dynamical systems. Recent experiments with 
confined two-dimensional electron gases (2DEG) revealed 
traces of chaos, nonlinear resonances, unstable periodic 
orbits or KAM-hierarchies in phase space 0, and also 
a variety of magnetic commensurability effects [USE! 
Among the latter, magnetic focusing is of particular 
importance and has found numerous applications such as 
the detection of composite fermions [5| or the separation 
of spin states . All these results were based on conduc- 
tance measurements. However, it is well known that the 
transition probabilities alone provide only an incomplete 
description of a scattering system. Important additional 
information is contained for example in the delay time, 
i.e., the time spent by the scattered particle in the in- 
teraction region. In a seminal work Wigner had pointed 
out the equivalence between the delay time and the en- 
ergy derivative of the phase of transition amplitudes • 
Ever since there was a lot of theoretical and experimen- 
tal activity devoted to understanding the distribution of 
this quantity in various physical contexts ||, including 
in particular also mesoscopic trans por t through systems 
with nonlinear dynamics |ft ll(lllllll2j . For semiconduc- 
tor nanostructures there exist pioneering measurements 
of picosecond delays for ballistic electrons in a magnetic 
focusing geometry Moreover, experimental access to 
the scattering phase (e. g., in quantum dots Q) might 
result in an electronic implementation of Wigner's rela- 
tion in the near future. 

Here we analyze the delay-time distribution for the 
transmission of electrons through a clean waveguide with 
a transversal magnetic field (Fig. QJi) . This model is ut- 
terly simple but not unrealistic in the above context and 
could be realized, e. g., as a finite constriction in a 2DEG. 
We find a rather surprising self-pulsing effect: A number 
of particles entering the interaction region at the same 
moment of time exit from the scatterer bunched together 
in packets which form a nearly periodic train. Specifi- 
cally we calculated numerically the distribution of trajec- 




FIG. 1: (a) A typical electron trajectory transmitting 
through a channel of length L under the influence of a perpen- 
dicular magnetic field, (b) Distribution of transmission delays 
for r = 4 and L = 200. The white dashed line is the average 
density, Eq. Q. The inset (c) magnifies the region around 
t* ~ 218.5 (Eq. IIOH 'I. Here the dashed vertical lines mark 
the delay times t n of orbits which are commensurate with the 
channel length (see text). The small insets demonstrate with 
two examples that the maxima t„ < t* (d) actually are log- 
arithmic singularities while this is not the case for t' n > t* 
(e). In these semilogarithmic plots the density is shown as a 
function of St — t — t n . The straight line in (d) has the slope 
given in Eq. 1251 . 



tory lengths from entry to exit and obtained a histogram 
which is composed of well-separated peaks on a smooth 
background (Fig.^J), c) (lif . 

Recently, self-pulsing was predicted in another physical 
context [ll|, and experimental confirmation came from 
microwave scattering in 2D resonators [J^- There the 
scattering echoes revealed a characteristic frequency of 
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the internal weakly chaotic dynamics. No similar expla- 
nation applies to our (integrable) model. Another strik- 
ing difference is the fact that in our case the delay-time 
density actually diverges at some of the peaks (Fig. |IJi). 

Our numerical results are reminiscent of transverse 
magnetic focusing where the conductance between 
two point contacts to a 2DEG oscillates as a function of 
the magnetic field (or some other parameter) and shows 
a regular pattern of singularities [17| . The fundamental 
dynamical origin of these singularities are caustics of the 
classical electron motion. They correspond to the sta- 
tionary points of a function that maps the incident angle 
to the exit location for an electron beam. Similarly we 
will explain the self-pulsing by stationary points of the 
delay time in the asymptotic phase space of scattering 
trajectories. To our knowledge, the effect of such sta- 
tionary points has never been analyzed before, not in any 
kind of scattering system. Therefore our results might be 
of interest beyond the study of electrons in a magnetic 
field, just as caustics in configuration space are relevant 
beyond magnetic focusing and in many different physical 
systems |l5j . We shall come back to this point at the end 
of the paper. 

We consider electrons entering the channel of Fig. 
from the left and with a given Fermi energy Ep = 
(m/2)v F . Otherwise the initial conditions < y < 1 and 
\tp\ < +-| are randomly drawn from the microcanonical 
distribution restricted to x — 0, p(y,<p) — \ cos ip. We 
use the channel width and the Fermi velocity as units 
of length and velocity, respectively (vf = b = 1). The 
cyclotron radius r = mvF/eB and the length L of the 
system are free parameters Any trajectory is com- 
posed of circular arcs which are traversed clockwise (by 
convention about the direction of the magnetic field) . For 
a point [x, y, ip) in phase space, the center of the current 
arc is the point (x c , y c ) — {x + r sin ip,y — r cos ip). Upon 
reflection from the channel walls y c is conserved. Thus 
the system is integrable. We restrict attention to y c < 0. 
For these trajectories the x-component of the velocity is 
positive at any moment of time, and they generate the 
bulk of the left-to-right transmission delay distribution 
(up to t ~ 1270 for the parameters of Fig.QJ. We change 
variables to (£,</>), where < <f> — arccos(— y c /r) < ir/2 
is the direction of the trajectory immediately after a re- 
flection from the lower wall. In terms of this angle, the 
average transport velocity of a trajectory in an infinite 
channel is v x = Ax(0)/At(</>). Here 



Ax{<p) = 2r(sin<£ — sin<£>) , 
At((j)) = 2r (0-$) 



(1) 
(2) 



are the horizontal and the total length of the trajectory 
segment between two reflections from the lower wall, re- 
spectively. If cos cf> < 1 — 1/r, the trajectory reaches 
the upper wall and is reflected there with the angle 
<f> = arccos(cos <j> + (see Fig. QJi). $ is set to zero for 



trajectories skipping along the lower wall only. £ is the 
longitudinal location of the first reflection from the lower 
wall. We make the convention — Ax ((/))/ 2 < £ < Ax{<p)/2 
such that negative values of £ correspond to positive tp 
and vice versa. Explicitly, the transformation between 
(<P, y) and (</>, £) is given by 

y = r(cos ip — cos <p) , (3) 
-2f = 2r(sin^-sin^) (f < 0) , (4) 
+2£ = 2r(sm0 + sin<p) (£ > 0) . (5) 

The probability density in the new variables is p{4>,£,) = 
\ sin </>. 

As a first approximation to the delay time we ignore 
the dependence on £ and set 



t(4>) = — = L 

v x Ax{(p) 

This results in the distribution 



p(t) = I d^d^ P (4,,os(t 

n0 t Ax((j> t ) 
2\T'U t )\ ' 



(6) 



(7) 



where <pt is the root of r(0) = t. For orbits which do not 
reach the upper wall we find 



r'W = L 

r 

P(t) = 



SIR <p t 



L sine 



k cos <p t 



(cos0 > 1-r- 1 ), (8) 
(L < t < t*) . (9) 



Clearly, the transmission delay is bounded from below by 
t = L. This corresponds to trajectories skipping along 
the lower wall with an infinitesimal reflection angle <j> — > 
such that their velocity is always directed along the 
channel. Close to the minimum we have according to 
Eq. © t - L(l + $/6) and from Eq. © P(t) - 3r(f> t /L. 
This explains the square-root like onset of the density, 
P(t) oc y/t — L, which is observed in Fig.^ 3 at t = 200. 
On the other hand the time 



t*(L,r) = L 



arccos(l — r x ) 
\/2r _1 — r -2 



(10) 



is the maximum delay for skipping orbits. Their con- 
tribution to P(t) abruptly drops to zero at this point 
(t* w 218.5 in Fig. QJj). Higher values t > t* corre- 
spond to trajectories which bounce off both, the lower 
and the upper wall of the waveguide. We omit here the 
lengthy expression replacing Eq. (JSJ) in this case. After 
substitution into Eq. J7J it yields the averaged delay-time 
distribution for t > t* which is displayed together with 
Eq. © by the dashed white line in Fig. \\}p. 

Now we come to the conspicuous oscillations in the 
delay-time density around its average value (Fig.^;). In 
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general the delay time t(6,^) depends on both variables, 
6 and £. A variation of £ corresponds to a longitudinal 
shift of the trajectory (Fig. [IJi) and results in [ltj 



d^r{6, £) = (cos <p) * — (cosy/) 



(11) 



where and ip' are the inclination angles of the trajectory 
on entry and exit, respectively. However, if the length 
of the system is an integer multiple of the length of a 
segment, 



L = nAx(6 n ) , 



(12) 



these angles are equal. Only in this case the dependence 
of the delay time on £ disappears and we have 



T(<pn,£) = r{(j> n ) = nAt(S n ) 



(13) 



For each n there can be at most two solutions to Eq. (|12|l . 
one corresponding to a skipping orbit and the other one 
not. We denote the corresponding delay times by t n < t* 
and t' n > t* , respectively. In Fig. QJ; these values are 
marked with dashed vertical lines, and we see that they 
coincide with the peaks in the delay-time distribution. It 
is tempting to explain this observation by the following 
simple (and misleading) argument: Typically, for a given 
6 the various values of £ contribute at different times 
to the density, while all of them add up in the same 
bin when 6 — (f> n . So it is clear that the count should 
have a peak there! However, this cannot explain why 
the peaks at t n are singularities of the density while the 
t' n are finite maxima. To highlight this fact numerically 
we have magnified in Fig. ^1 the density in the vicinity 
of £40 and observe that it grows as P(t) oc ln(f — t n ). In 
contrast the corresponding plot for t' i0 (Fig. QJ) saturates 
to a finite value. 

In order to understand the nature of the maxima we 
construct the function r(0, £) explicitly and represent the 
delay-time density as 



P(t) = Jdd fd£p(<l>,Z)8(t-T{tl>,t)) 



(14) 



<j>=Ht,£) 



To this end we decompose a trajectory into a number 
n of complete arcs and two additional terms for entry 
and exit. Formally we define n(6) as the best integer 
approximation to L/Ax, e. g. n = 3 in Fig. The 
mismatch to exact commensurability will be denoted by 

8L(6) = nAx(6) — L . (15) 

According to Fig. ^ the delay time is given by 

t(0,0 = nAt(0)+a(0,O-*(0,O (16) 

with C = £ + SL{6). In Eq. (JTSJl, s(<f>,£) is the length 
of the incomplete arc segment at the beginning of the 



trajectory. For £ in Fig. ^1 this is 

s((j>, f) = r{6 + if) {ip < 0, < £ < Ax(6)/2) 

= r[6 + arcsin(£/r — sin (f>)] (17) 

(the second line follows from Eq. JSJl). However, £' can 
take values in the interval [— Ax (6), Ax(6)] since we have 
\8L(6)\ < £f and |£| < As. Therefore we extend the 
domain of s(6, £) by the definitions 



s(0, Ax(6) 



-0 





Ai(0)- s (0,O. 



(18) 
(19) 



Note that s(<f>,£) is continuous at s((j), 0) = and 
s(6,Ax(6)/2) = At{6)/2 but not analytic at £ = 0. 
Eqs. JTjj-lO represent the delay time for < 6 < f 
and arbitrary values — A^ < £ < +A^. Fig. |3 displays 
this function in the regions (a) t 4 i < i < t 39 and (b) 
^41^;^^ *39- All level sets of constant delay time 
are one-dimensional curves in the two-dimensional phase 
space. In this respect there is nothing special about 
t = t n . Rather we observe in Fig. 0: c^t(0„,O) = 0. 
This is not immediately obvious from geometrical con- 
siderations but it can be confirmed analytically. We con- 
clude that (</>,£) = (6 n ,0) are stationary points of the 
delay time, and this is what really singles out the values 
t = t n . No stationary points exist for t > t* (Fig. EJl), 
although d^r is very small close to £ = ±Ax(6)/2. 
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FIG. 2: Lines of constant delay time t(£, (j>) =t are shown in 
the vicinity of (a) t^o and (b) t' i0 . The function d^r vanishes 
at (j> = (j> n and £ = (c) while no stationary point exists for 

4> = Cp'n (d). 

In the vicinity of the stationary points we can expand 
the delay time to second order, 



St = Id^rS^ + d^r^Scf, + l -dl C e , (20) 



solve for 56 = 



and substitute the result into 



(21) 



This gives d^r as a function of £ and St = t—t n and can be 
used in Eq. (|14f) to calculate the divergent contribution to 
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the delay-time density. Unfortunately this procedure is 
more complicated than it may seem on first sight. Due to 
Eq. (|18|l we must distinguish for St, 5L > three regions 
(i) £ > 0, (ii) £ < £ < and (iii) £ < £ ( an d again three 
regions for St, SL < 0) . Here £ (<^) is the point where the 
third term in Eq. changes sign. It is given implicitly 
by the condition £o = —SL((f> n + S(f>) where S<j)({;o,St) is 
the root of Eq. 1121 )[l . Solving for £o we find 



£ = - (sin 4> n ) - l ^28t (cos 3 cf> n ) / (2n - 1 ) 



(22) 



Note that we did not drop the third term in Eq. (|20|l . 
In fact we have 9|^r 7^ in the second region although 
r(0„, £) is constant as a function of £. This is no contra- 
diction since, according to Eq. (|22p. this region shrinks 
to a point for 4> = 4> n (and hence is not visible in Fig. EJi) . 
With Eqs. (|20|I - H22() it can be shown that the contribution 
of region (ii) approaches a constant as St — > 0. Hence it 
is irrelevant for our purpose. In the regions (i), (iii) we 



use d, 



and find 



d+r&St) = J(d, 



r)H 2 + 2Std^r 



(23) 



1^ 



In + \d+T\) + const. (24) 



t n )\ (t->tn). (25) 



This prediction for the strength of the divergence has 
been confirmed numerically. For example both, a fit to 
the data in Fig.[TJi and Eq. yield a prefactor 0.00381 
for n = 40. 

The rather complicated analytical structure in the 
vicinity of the stationary points of the delay time, and 
also the self-pulsing due to their regular arrangement, are 
specific to our model. We would like to stress, however, 
that these features are not at all necessary for the appear- 
ance of singularities in the delay-time density. Generic 
stationary points for 2D scattering problems are extrema 
and saddle points. It is easy to see that saddle points 
lead again to a logarithmic divergence, while the den- 
sity is finite at isolated maxima and minima. Although 
a logarithm is a rather weak singularity, it should lead 
to a prominent maximum in a histogram with finite res- 
olution such as Fig.^J:. This might be important, e. g., 
for inverse scattering problems. Consider for example a 
situation where it is difficult to control the precise initial 
conditions of test particles while their time delay can be 



measured with high precision. Then saddle points of the 
delay time will be among the dynamical features which 
are directly accessible from experimental data. There- 
fore it seems worthwhile to study them in more generic 
situations and independent of self-pulsing effects. 
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This is not the case in the absence of a magnetic field (r = 
00). There the distribution is monotonically decreasing, 
Pit) = L- x {L/tf{\ - [L/tf)- 1/2 for t > L. 
In reality these singularities are blurred by the finite 
contact width or by quantum effects 0. An analogous 
smoothing of the delay-time singularities must be ex- 
pected, but this will not be considered here. 
Our results are not sensitive to the precise values of L and 
r as long as the magnetic field has a moderate strength 
r > 1. For numerical calculations we choose r = 4 and 
L = 200. 

We use a notation where, e. g., d^r stands for dr/d^. All 
second-order derivatives are evaluated at the stationary 
point (0,£) = (0 n ,O). 



